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Abstract. We discuss the exclusive radiative B-meson decays B — > K*y and B — > py in a model- 
independent manner. The analysis is based on the heavy-quark limit of QCD. This allows a factoriza- 
tion of perturbatively calculable contributions to the B — > Vy matrix elements from non-perturbative 
form factors and universal light-cone distribution amplitudes. These results allow us to compute ex- 
clusive b — ► s(d)y decays systematically beyond the leading logarithmic approximation. We present 
results for these decays complete to next-to-leading order in QCD and to leading order in the heavy- 
quark limit. Phenomenological implications for various observables of interest are discussed, in- 
cluding direct CP violation and isospin breaking effects. 

INTRODUCTION 

The radiative transitions b — > s(d)y are among the most valuable probes of flavour 
physics. Although they are rare decays the Cabibbo-favoured b — > sy modes are ex- 
perimentally accessible already at present. The inclusive branching fraction has been 
measured to be 

B(B -> X s y) = (3.23 ±0.42)- 10~ 4 (1) 

combining the results of [2, 3, 4]. The branching ratios for the exclusive channels have 
been determined by CLEO [5], and more recently also by BaBar [6] and BELLE [7] to 
give the averaged values: 

B(B° -> K*°y) = (4.51 ±0.54) • 10~ 5 (2) 
B(B + ^K* + y) = (3.86 ±0.59)- 10~ 5 . (3) 

On the theoretical side, the flavour-changing neutral current (FCNC) reactions b — > s(d)y 
are characterized by their high sensitivity to New Physics and by the particularly large 
impact of short-distance QCD corrections [8, 9, 10, 11]. Considerable efforts have 
therefore been devoted to achieve a full calculation of the inclusive decay b — > at next- 
to-leading order (NLO) in renormalization group (RG) improved perturbation theory 
[12, 13, 14] (see [15] for recent reviews). 
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Whereas the inclusive mode can be computed perturbatively, using the fact that the 
b-quark mass is large and employing the heavy-quark expansion, the treatment of the 
exclusive channel B — > K*y is in general more complicated. In this case bound state 
effects are essential and need to be described by nonperturbative hadronic quantities 
(form factors). The basic mechanisms at next-to-leading order were already discussed 
previously for the B — > Vy amplitudes [16]. However, hadronic models were used to 
evaluate the various contributions, which did not allow a clear separation of short- 
and long-distance dynamics and a clean distinction of model-dependent and model- 
independent features. 

In this talk we present the results of [1] where a systematic analysis of the exclusive 
radiative decays B — > Vy (V = K*, p) in QCD, based on the heavy quark limit mt, 3> 
Agco* was performed. We quote factorization formulas for the evaluation of the relevant 
hadronic matrix elements of local operators in the weak Hamiltonian. A similar subject 
was treated in [17, 18]. Factorization holds in QCD to leading power in the heavy quark 
limit. This result relies on arguments similar to those used previously to demonstrate 
QCD factorization for hadronic two-body modes of the type B — > nn [19, 20]. 

This framework allows us to separate perturbatively calculable contributions from 
the nonperturbative form factors and universal meson light-cone distribution amplitudes 
(LCDA) in a systematic way. This includes the treatment of loop effects from light 
quarks, in particular up and charm. Such loop effects are straightforwardly included 
for the inclusive decays b — > s(d)y. For the exclusive modes, however, the effects from 
virtual charm and up quarks have so far been considered to be uncalculable "long- 
distance" contributions and have never been treated in a model independent fashion. 

Finally, power counting in Agco/"^ implies a hierarchy among the possible mech- 
anisms for B — > Vy transitions. This allows us to identify leading and subleading con- 
tributions. For example, weak annihilation contributes only at subleading power in the 
heavy quark limit. 

Within this approach, higher order QCD corrections can be consistently taken into 
account. We give the B — > Vy decay amplitudes at next-to-leading order (NLO). Fur- 
thermore numerical values for CP-asymmetries and isospin-violating ratios are given. 
After including NLO corrections the largest uncertainties still come from the B — > V 
form factors, which are at present known only with limited precision (~ ±15%), mostly 
from QCD sum rule calculations [21]. The situation should improve in the future with 
the help of both lattice QCD [22] and analytical methods based on the heavy-quark and 
large-energy limits [23, 24, 25]. 



BASIC FORMULAS 



The effective Hamiltonian for b — > sy transitions reads 



V ^ p=u,c 



CiQ p l +C 2 Q p 2 + £ GQi 



i=3,...,8 



(4) 



FIGURE 1. Contribution of the magnetic penguin operator Qi described by B — > V form factors. All 
possible gluon exchanges between the quark lines are included in the form factors and have not been 
drawn explicitly. 



(s) 

where Xp = V* s V p b- The relevant operators are given by 

Q\ = (sp) V -A(pb)v-A (5) 

Q P 2 = (SiPj)v-A(Pjbi)v-A (6) 

Qi = ^mbSiC^il+ys^iF^ (7) 
Qs = ■^m b s i & JV (l+y 5 )T i a j bjG^ (8) 

Note that the numbering of Q\ 2 is reversed with respect to the convention of [26]. We 
neglect the contribution from the QCD penguin operators Q^. ..6, which enter at 0(a s ) 
and are further suppressed by very small Wilson coefficients. The effective Hamiltonian 
for b — > dyis obtained from (4-8) by the replacement s — > d. 

The most difficult step in computing the decay amplitudes is the evaluation of the 
hadronic matrix elements of the operators in (4). A systematic treatment can be given in 
the heavy-quark limit. In this case the following factorization formula is valid 



(Vy(e)\Qi\B) = 



F B^V (0) T l + 1 1 ^ dv T I1 £ y) ^ ^ ^ (y) 

Jo 



(9) 



where £ is the photon polarization 4- vector. Here F B ^ V is a B — > V transition form factor, 
and <J»b, <Ity are leading twist light-cone distribution amplitudes of the B meson and the 
vector meson V, respectively. These quantities are universal, nonperturbative objects. 
They describe the long-distance dynamics of the matrix elements, which is factorized 
from the perturbative, short-distance interactions expressed in the hard- scattering kernels 
T/ and T/ 1 . The QCD factorization formula (9) holds up to corrections of relative order 

In the leading logarithmic approximation (LO) and to leading power in the heavy- 
quark limit, Qi gives the only contribution to the amplitude of B — > Vy and the factor- 
ization formula (9) is trivial. The matrix element is simply expressed in terms of the 
standard form factor, T 7 7 is a purely kinematical function and the spectator term T 7 /7 is 
absent. An illustration is given in Fig. 1. The matrix element reads 



(V(k,r\)y(q,E)\Q 7 \B) = -^m h c v F v \E^ x %r\ v k x q p + i(e-r]k-q-E-kr]-q) 



(10) 





FIGURE 2. o(0Cs) contribution to the hard-scattering kernels T( from four-quark operators Q\. The 
crosses indicate the places where the emitted photon can be attached. 
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FIGURE 3. O (a. s ) contribution to the hard-scattering kernels r 8 7 from chromomagnetic penguin opera- 
tor g 8 . 



where c v = 1 for V = K*, p and c v = 1 /V2 for V = p°. The B — > V form factor Fy is 
evaluated at momentum transfer g 2 = 0. Our phase conventions coincide with those of 
[21,27]. 

The matrix elements of Q\ and Qg start contributing at 0(a s ). In this case the 
factorization formula becomes nontrivial. The diagrams for the hard- scattering kernels 
T/ are shown in Fig. 2 for Q i and in Fig. 3 for Q$. These diagrams were computed 
in [13] to get the virtual corrections to the inclusive matrix elements of Q\ and Qg. 
In our case they determine the kernels T\ and T%. As required for the consistency of the 
factorization formula these corrections must be dominated by hard scales of order m b and 
hence must be infrared finite. This is indeed the case. Re-interpreted as the perturbative 
hard-scattering kernels for the exclusive process, the results from [13] imply 

(Qi,*)' = (to^Gifi (11) 

where Cf = (N 2 — 1)/ (2N), with Af = 3 the number of colours, and 

Gi(s c ) = -^ln-^-+ gl (s c ) (12) 
II nit, 

G 8 = ^ln^-+^ 8 (13) 
3 mb 

2 

The finite part gi(s c ) of the Q i -contribution depends via s c = ^§ on the mass of the 

quark running in the loop in Fig. 2. 

We now turn to the mechanism where the spectator participates in the hard scattering. 
The non-vanishing contributions to T t n are shown in Fig. 4. To find the correction for 




FIGURE 5. The hard-scattering kernel h{v,s c ) as a function of v. 



(<2i) we compute the first diagram in Fig. 4 We obtain 

{Ql f = { Q 7) C ^pL Hl ( Sc ) (14) 

with 

The first negative moment of the light-cone distribution amplitude <5bi(£) can be pa- 
rametrized by mg/^e where = C(Agco)- The hard-scattering kernel h(v,s c ) is dis- 
played in Fig. 5. It is real for v <As and develops an imaginary part for the light-cone 
momentum fraction of the anti-quark in the vector meson v > As. 

The correction to (Qg) from the hard spectator interaction comes from the second 
diagram in Fig. 4. One finds 

{Q ,f = {Qi) ^M£Lh 8 (16) 



where 

i B Jo 



H ^M^f d Mf dv ^l (17) 



3N F v m 2 B Jo ^ 



Here the LCDA (v) is of leading power for a transversly polarized vector meson. The 
last convolution integral in (17) can be performed explicitely and leads to a combination 
of Gegenbauer moments of [28, 21]. 

There are further mechanisms that can in principle contribute to B — > Vy decays. One 
possibility is weak annihilation. In this case the leading-power projection onto the meson 
V vanishes because the trace over an odd number of Dirac matrices is zero. A non- 
vanishing result arises from the subleading-power projections onto <J>y and 

Despite its power suppression, the dominant annihilation amplitude can be computed 
within QCD factorization. This is because the colour-transparency argument applies to 
the emitted, highly energetic vector meson in the heavy-quark limit [20, 29]. 

Since weak annihilation is a power correction, we will content ourselves with the 
lowest order result (o(cc^)) for our estimates below. In particular, we shall include 
the annihilation effects from operators <2i,2 to estimate i so spin-breaking corrections in 
B — > py decays. The reason for including this class of power corrections is that they come 
with a numerical enhancement from the large Wilson coefficients (C\ ^3\Cq\) and 
are not CKM suppressed. Instead, a CKM suppression of annihilation effects occurs for 
B — > K*y and these contributions are thus very small in this case. 

RESULTS 

Finally, we can combine these results and write, adding the up- and charm- quark 
contribution, 

(Vy\Q 7 \B) (18) 

where, at NLO, the factorization coefficients of (Vy) are given as 

a p 7 (Vy) = C 7 + ^^(C^G^+Cs^Gs) 

+ asi ^ CF (CMH^+CMHs) (19) 

Here the NLO expression for C-j [14] has to be used, while the leading order values are 
sufficient for C\ and Cg. Numerically we obtain for central values of all input parameters, 
at fx = nib, and displaying separately the size of the various correction terms: 

a c 7 (K*y) = -0.3221 +0.0113 -0.0820-0.0147? -0.0144-0.0109/ 
dj° AC^ L0 r/ 8 -contribution r/ 7 8 -contribution 

= -0.4072 -0.0256/. (20) 

We note a sizable enhancement of the leading order value, dominated by the r 7 -type 
correction. A complex phase is generated at NLO, where the T^-corrections and the 
hard- spectator interactions (T n ) yield comparable effects. 

The net enhancement of a-j at NLO leads to a corresponding enhancement of the 
branching ratios, for fixed value of the form factor. This is illustrated in Fig. 6, where we 
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FIGURE 6. Dependence of the branching fractions B(B° — > K*°y) and B(B~ — > p~y) on the renor- 
malization scale /j. The dotted line shows the LO, the dash-dotted line the NLO result including type-I 
corrections only and the solid line shows the complete NLO result. 



show the residual scale dependence for B(B —> K*°y) and B(B — > p~y) at leading and 
next-to-leading order. As expected, the inclusion of the hard-vertex corrections (T 1 ) re- 
duces the scale dependence coming from the Wilson coefficients. The scale dependence 
of the complete NLO result is "deteriorated" because the type-II contributions appear at 
O (a s ) for the first time and therefore introduce a completely new scale dependence. 

For the decay B — > pyboth sectors of the effective Hamiltonian have the same order of 
magnitude. The amplitude for the CP-conjugated mode B — > py is obtained by replacing 

)Sp^ — > Xp^*. We may then consider the CP asymmetry 

^ (py) " r( B ^p Y)+ r(^py) <21) 

It is substantial for the py modes and much less dependent on the form factors. 
Here the largest theoretical uncertainty comes from the scale dependence. This is to 
be expected because the direct CP asymmetry is proportional to the perturbative strong 
phase difference, which arises at 0(a s ) for the first time. Unknown power corrections 
could have some impact on the prediction. 

A further interesting observable is the charge averaged isospin breaking ratio 

Afnri r(* + -> P + Y) | r(*~ -> p-y) 



Within our approximations, isospin breaking is generated by weak annihilation. Isospin 
breaking was already discussed in [30], partially including NLO corrections. 
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FIGURE 7. Left: The CP asymmetry ^cp(PY) as a function of the CKM angle Y for three values of the 
renormalization scale /a = nib/2, nib an d 2m£. Right: The isospin-breaking asymmetry A(py) as a function 
of the CKM angle Y at leading and next-to-leading order. 



B — > pyalso depends sensitively on fundamental CKM parameters, such as \V u b/V c t,\ 
and y, and can thus in principle serve to constrain the latter quantities once measurements 
become available. This is further illustrated in Fig. 7, where the dependence on y is 
shown for Acp(py) and A(py), respectively. We remark that our sign of A(py) differs 
from the one found in [30]. 

Another application of our results concerns an estimate of U-spin breaking effects in 
B — > Vy decays [31, 32, 33]. Let us state here only that U-spin breaking effects can be 
sizeable. For further details we refer to [1]. 



In this talk we have presented a systematic and model-independent framework for the 
exclusive radiative decays B — > Vy based on the heavy-quark limit. This allowed us to 
compute the decay amplitudes for these modes consistently at next-to-leading order in 
QCD. 

An important conceptual aspect of this analysis is the interpretation of loop contri- 
butions with charm and up quarks, which come from leading operators in the effective 
weak Hamiltonian. We have argued that these effects are calculable in terms of pertur- 
bative hard-scattering functions and universal meson light-cone distribution amplitudes. 
They are 0(cc s ) corrections, but are leading power contributions in the framework of 
QCD factorization. This picture is in contrast to the common notion that considers charm 
and up-quark loop effects as generic, uncalculable long-distance contributions. Non- 
factorizable long-distance corrections may still exist, but they are power-suppressed. 

Another important feature of the NLO calculation are the strong interaction phases, 
which are calculable at leading power. They play a crucial role for CP violating observ- 
ables. We have seen that weak-annihilation amplitudes are power-suppressed, but can 
be numerically important for B — > py because they enter with large coefficients. These 



CONCLUSIONS 



effects also turn out to be calculable and were included in our phenomenological discus- 
sion at leading order in QCD. 

Finally, a numerical analysis leads to our predictions for the branching ratios B(B° — > 
K*°y) ~ 7.09 • 1(T 5 and B(B -> p"y) ~ 1.58 • 1(T 6 , and for the CP asymmetries 
Acp(K*y) ~ —0.5% and ^cp(PY) ~ 10% for central values of the input parameters. 
The uncertainties for the branching ratios are of 0(30. . .35%) where the form factors 
Fk* and F p taken from [21] clearly dominate this uncertainty. This situation however 
can be systematically improved. In particular, our approach allows for a consistent 
perturbative matching of the nonperturbative form factor to the short-distance part of 
the amplitude. Power corrections are not expected to be unreasonably large because 
with the annihilation contribution we have included the probably numerically largest of 
these corrections. Furthermore there are no chirally enhanced power corrections in the 
vector meson wave functions. 
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